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We investigate the QCD Dirac operator with gauge configurations given by a liquid of instantons in the 
region of temperatures about the chiral phase transition. Both the quenched and unquenched cases are examined 
in detail. We present evidence of a localization transition in the low lying modes of the Dirac operator around 
the same temperature as the chiral phase transition. Thus both level statistics and eigenvectors of the QCD 
Dirac operator at the chiral phase transition have similar properties than those of a disordered conductor at the 
metal-insulator transition. This strongly suggests the phenomenon of Anderson localization (localization by 
destructive quantum interference) is the leading physical mechanism in the restoration of the chiral symmetry. 
Finally we argue that our findings are not in principle restricted to the ILM approximation and may also be 
found in lattice simulations. 
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The spontaneous breaking of the approximate chiral sym- 
metry (S^SB) and its eventual restoration at finite temperature 
is one of the most important features of the strong interactions. 
The order parameter associated to this symmetry breaking is 
the chiral condensate, (vj/\|/) . By definition, the condensate is 
nothing but a quark loop in momentum space. Thus it should 
naively vanish as the quark mass goes to zero. In nature the 
lightest quarks are not massless so a nonzero condensate is 
expected even in a free theory. However the small quark bare 
mass can only account for a very small percentage of the chi- 
ral condensate, the rest has its origin in the strong non pertur- 
bative color interactions of QCD. QCD models whose gauge 
configurations are given by an interacting liquid of instantons 
JH] provide an adequate theoretical framework to understand 
S^SB. Thus, based on the semiclassical picture of a QCD vac- 
uum dominated by instantons, it has been suggested yfl that 
the S^SB in QCD and the phenomenon of conductivity may 
have similar physical origins. Conductivity in a disordered 
sample is produced by electrons that although initially bound 
to an impurity may become delocalized by orbital overlap- 
ping with nearby impurities. Similarly, in the QCD vacuum, 
the zero modes of the Dirac operator though initially bound 
to an instanton may get delocalized due to the strong overlap 
with other instantons. As a consequence, the chiral conden- 
sate becomes nonzero and chiral symmetry is spontaneously 
broken. In this paper we provide evidence that these analo- 
gies can be extended to describe the chiral phase transition 
in QCD at finite temperature. Specifically, in the context of 
an Instanton Liquid Model (ILM), we show that around the 
temperature that the chiral phase transition occurs the low ly- 
ing eigenmodes of the QCD Dirac operator undergo an An- 
derson transition (AT), namely, a localization-delocalization 
transition, characterized by multifractal eigenstates and criti- 
cal statistics |4J]. 

For the sake of clarity we provide a brief summary of the 
properties of a disordered system at an AT. We recall that in 
three and higher dimensions there exists a mobility edge sep- 



arating localized from delocalized states. In the delocalized 
region, eigenfunctions are extended through the sample and 
level statistics are described by random matrix theory (RMT). 
In the opposite limit, eigenfunctions are exponentially local- 
ized and the spectral correlations are described by Poisson 
statistics. Around the mobility edge eigenstates are multifrac- 
tal |0] and level statistics are described by critical statistics 
10]. Typical features of critical statistics include scale invari- 
ant spectrum 111 ill , level repulsion and asymptotically linear 
number variance. 

Thus the level spacing distribution, P(s) (the probability of 
having two eigenvalues at a distance s) goes to zero as the 
spacing s separation does, P(s) — > s — > 0. The number vari- 
ance T?(l) = {(N e - (Ne)) 2 } ~ %£ for I > 1 (Ne is the number 
of eigenvalues in an interval of length £) is asymptotically lin- 
ear, as for an insulator (% = 1), but with a slope % < 1 (~ 0.27 
for a 3D Anderson model [13]). Below we show that at the 
chiral restoration temperature both eigenvalues and eigenvec- 
tors of the QCD Dirac operator have similar properties. 



I. INSTANTON LIQUID MODEL 

In this section we highlight the main features of the QCD 
vacuum model based on instantons to be used in this paper. 

Instantons [3] are classical solutions of the Yang-Mills 
equations of motion which minimize the action in Euclidean 
space. They are believed to be the leading semiclassical con- 
tribution to the bosonic part of the QCD path integral. How- 
ever the construction of a consistent QCD vacuum based on 
instantons faces serious technical difficulties. Exact multi- 
instanton solutions are hard to obtain since the Yang-Mills 
equations of motion for QCD are nonlinear and therefore a 
superposition of single instanton contributions is not itself 
a solution. Additionally quantum corrections may spoil the 
semiclassical picture implicitly assumed of a QCD vacuum 
composed of instantons well separated and weakly interact- 
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ing. These problems have been overcome either by invok- 
ing variational principles Q5|] or by phenomenologically fixing 
certain parameters of the instanton ensemble. The latter case, 
usually referred to as the instanton liquid model (ILM) 1 1 ] 
(for a modern review see II 1510 . yields accurate estimates of 
vacuum condensates and hadronic correlation functions ifToh 
for most light hadrons just by setting the mean distance be- 
tween instantons to be R w 1 fm (which corresponds to a den- 
sity N/V ~ 1 fm~ 4 ) and the mean instanton size psj 1 /3 fm. 
Lattice simulations have also supported the picture of a QCD 
vacuum dominated bv instantons I U4I 1 . 
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FIG. 1: Level spacing distribution P(s) in the bulk for the quenched 
ILM at T = 200 MeV for different volumes. The inset shows logP(s) 
for the tail of the data. The best fit (solid line) corresponds to a slope 
-1.64. 
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FIG. 2: Spectral rigidity in the bulk for the quenched ILM at T = 200 
MeV for different volumes. The result has very little size dependence 
and agrees well with the prediction of critical statistics (solid line). 



A. Technical details of the numerical simulation 

The ILM partition function for Nf quark flavors with 
masses m% is given by 
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where the integral is over the positions, sizes and orienta- 
tions of the instantons and Sym is the Yang-Mills action. The 
fermion determinant is evaluated in the space of the fermionic 
zero modes of the instantons. For further discussion of this 
partition function we refer to 01511 . We just mention that we 
use the standard phenomenological value of the instanton den- 
sity N/V = lfm~ 4 . For the sake of simplicity we have kept 
this density fixed as the temperature is increased. We justify 
this approximation based on the fact that even for tempera- 
tures close to the chiral phase transition the instanton density 
is still sizable (around 0.6 fm~ 4 lfl5lP . Indeed in previous sim- 
ulations the drop in the density was ruled out as the physical 
mechanism leading to the chiral phase transition jl5ll . 

All units in the ILM are typically given in terms of the 
QCD scale parameter A which we have simply set to 200 
MeV. This choice is sufficient for our purposes since we are 
not concerned with making quantitative predictions about the 
position of the QCD chiral phase transition, but instead are 
interested in more general features of the transition that we 
don't expect to be very sensitive to a fine tuning of the param- 
eters of the ILM. We therefore use the ILM as a qualitative 
model for QCD at temperatures around the chiral phase tran- 
sition. We also stress that it provides a reasonable description 
of S%SB and many hadronic correlation functions both at zero 
and at finite temperature II 1 C XI . Furthermore the ILM has mod- 
est computational requirements and allows us to go to fairly 
large volumes and get good statistics. 

In this paper we present results for both the quenched 
Nf = and the unquenched Nf = 2 cases (Nf = 2 with m u = 
ma = 0). The partition function, equation (Q]), is evaluated us- 
ing a standard Metropolis algorithm. In the quenched case we 
performed between 2000 to 10000 measurement sweeps for 
each set of parameters after allowing around 1000 sweeps for 
thermalization. Results are presented for ensembles of up to 
N = 5000 instanton and anti-instantons. In the unquenched 
case we were only able to investigate ensembles of up to 
N = 700 instantons and anti-instantons. We also discarded the 
first 1000 sweeps in each simulation and did between 1000 to 
5000 sweeps per ensemble. Thermodynamic properties of the 
model have already been addressed in lfl5tl and will not be 
discussed here. We will focus mainly on observables such as 
level statistics and eigenvector scaling properties of the QCD 
Dirac operator which are especially well suited for studying 
Anderson localization effects 111 811 . All the spectral correlators 
are calculated from the unfolded spectrum. This procedure 
scales the eigenvalues so that the spectral density on a spec- 
tral window comprising several level spacings is unity but it 
does not remove the small fluctuations about the mean density 
that provide fundamental information about the system. 

Finally we recall that in the ILM the effective disorder pa- 



2 



Brazilian Journal of Physics, vol. 37, 246, (2007) 



rameter is the temperature since we have fixed the density of 
instantons. One of the most challenging tasks of the numeri- 
cal calculation is to accurately locate the temperature at which 
the AT occurs. In all cases in the paper this temperature has 
been estimated by the finite size scaling method introduced 
in [11] in the context of disordered systems. In essence this 
method consists of computing a spectral correlator, such as the 
level spacing distribution P(s) or the number variance T. 2 (£), 
for different sizes and then finding the temperature at which it 
becomes size independent. 



II. LOCALIZATION IN THE INSTANTON LIQUID 
MODEL: THE QUENCHED CASE 

In this section we present results for the quenched ILM 
at nonzero temperature. Without dynamical fermions there 
is not technically a chiral phase transition since there is no 
chiral symmetry. However one can in principle study the 
"quenched" quark condensate in a purely gluonic background 
through the Banks-Casher relation 11611 . However, in the 
quenched ILM, the spectral density seems to diverge close to 
the origin 101 (similarly to recent lattice results with overlap 
fermions J21) thus suggesting a likewise divergent "quenched" 
quark condensate. In view of these facts, the transition to lo- 
calization (see below) we have observed in the quenched QCD 
Dirac operator at finite temperature both close to origin and in 
the bulk of the spectrum can't be linked directly to either of 
these phenomena. We still consider the quenched ILM results 
to be of interest as we have found a very clear example of a 
mobility edge with features strikingly similar to those of a 3D 
disordered system at the AT. 

First we examine the properties of eigenstates in the bulk 
of the spectrum. This region is not directly related to the chi- 
ral transition however it is a very clear example of a mobility 
edge in the ILM. Moreover since it affects a larger range of 
eigenmodes we can get better statistics. By using the finite 
scaling method introduced in 111 111 we have observed a mobil- 
ity edge separating localized from delocalized eigenstates in 
the range T ~ 150 - 250 MeV. As the temperature decreases 
the location of the mobility edge moves to the end of the spec- 
trum. For T < 170 MeV the results are less reliable since the 
mobility edge is located almost in the end of the spectrum 
where truncation effects are larger. Unlike the unquenched 
case, where the first temperature to study is in principle dic- 
tated by the restoration temperature, we do not have a special 
reason to choose a specific temperature to investigate the AT. 

Following the literature in disordered systems we have in- 
vestigated the temperature T ~ 200 MeV such that the mobil- 
ity edge is located around the center of the spectrum, namely, 
if the spectrum spans a spectral window [0, e/], the center cor- 
responds with a region around ~ ef/2. 

The level statistics in this spectral region have all the sig- 
natures of an AT including a spectrum that is scale invari- 
ant to high degree. As observed in figures Q] and [2] both 
short range statistics such as the level spacing distribution, 
P(s), and long range correlators such as the spectral rigidity 
A 3 (£) = if fo(£ 3 -2l 2 x + x 3 )Y?(x)dx do not depend on the 



system size for volumes (number of instantons) ranging from 
N = 500 to N = 5000. Level repulsion is still present as for 
a metal but the spectral rigidity (see figure |2} is asymptoti- 
cally linear with a slope that corresponds to % ~ 0.29 ±0.02 
in fair agreement with the value for a 3D disordered system at 
an AT of % ~ 0.27 ± 0.02 [ 1 3] . Furthermore the exact form of 
the spectral rigidity follows closely the prediction of critical 
statistics, equation (31) in [ 17] with h, a free parameter, set to 
a value of 0.62. 

We now investigate the multifractality of the eigenstates 
by looking at the scaling of P q = J d d r\y(r)\ 2q ^ L~ D i( q ~ l \ 
where D q is a set of different exponents describing the tran- 
sition and L °< N 1 ^ is the spatial system size. The second 
moment P2 is usually referred to as the inverse participation 
ratio (IPR). 

In order not to consider eigenstates outside the critical re- 
gion we have taken only 10% of the eigenvectors around the 
center of the spectrum for T — 200 MeV. By fitting P2 (also 
known as the inverse participation ratio or IPR) for different 
volumes, we have obtained a value of D2 ~ 1 .5 ± 0. 1 in agree- 
ment with that of a 3D disordered system at the AT 01811 . 

Our spectral and eigenfunction analysis shows that at T = 
200 MeV there is mobility edge around the central part of the 
spectrum in the quenched ILM. We have also observed simi- 
lar critical features in the region close to the origin at tempera- 
tures around 1 10-140 MeV. Here we expect the appearance of 
the mobility edge to be related to S^SB, However the analysis 
is complicated by the accumulation of very small eigenval- 
ues typical of the quenched approximation in QCD. Addition- 
ally numerical results are less conclusive due to poor statistics. 
Despite these limitations we have also observed critical statis- 
tics and multifractal eigenstates with D2 ~ 0.8 ±0.1 (more 
sparse than in the bulk). 
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FIG. 3: Chiral condensate and inverse participation ratio (IPR) of the 
lowest eigenmode for the ILM with two massless quarks. 
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FIG. 4: Level spacing distribution P(s) at the origin for the un- 
quenched ILM at T — 115 MeV for different system sizes. 



III. CHIRAL RESTORATION AND ANDERSON 
TRANSITION: THE UNQUENCHED CASE 

We now study eigenvalues and eigenvectors of the QCD 
Dirac operator in an ensemble with two massless quark fla- 
vors. In this case the chiral condensate associated with S^SB 
is finite. Around the critical temperature, the condensate ap- 
proaches zero signaling the chiral phase transition. Since we 
have two massless flavors we expect to find a second order 
chiral phase transition. 

We focus our attention on the spectral region close to ori- 
gin since our aim is to relate the chiral phase transition with 
a transition to localization (AT) of the QCD Dirac operator 
eigenmodes. By looking at the chiral condensate versus tem- 
perature for a range of system sizes (figure O we see that the 
condensate approaches zero around temperatures of 1 15 - 120 
MeV in a manner consistent with a second order phase tran- 
sition. At the same time the condensate is falling we see the 
IPR of the lowest eigenmode begin to rise signaling a transi- 
tion to localization. This finding strongly suggests that both 
phenomena are intimately related as was first conjectured in 
Ref.fl. 

Next we examine the nature of the localization transition. 
First we locate, by using the finite size scaling method Hill , 
the AT for the lowest lying eigenmodes at T ~ 115-1 20 MeV. 
At this temperature we see that eigenvalue statistics such as 
the level spacing distribution (figure|4|i are nearly scale invari- 
ant, there is level repulsion but the tail of P(s) is exponential 
as at the AT. We also looked at the scaling of eigenmode mo- 
ments P q in the lowest eigenmode at T = 120 MeV for system 
sizes ranging from I? =63-252. This yielded a set of fractal 
dimensions D 2 = 1.3 ± 0.2, D 3 = 0.9 ± 0.2, D 4 = 0.7 ± 0.3 
and D5 = 0.7 ± 0.3 Although higher volumes would be de- 
sirable to fully confirm the multifractal scaling of the eigen- 
states, it is encouraging that the IPR of comparatively similar 
temperature (T ~ 100 MeV) is comparable to the one for a 



fully metallic sample. We also mention that we observed a 
mobility edge on the bulk of the spectrum at a slightly higher 
temperature with properties similar to the quenched case. 

At this point a natural question to address is whether the 
results here reported concerning the Anderson transition in the 
ILM may also be present in more realistic models of the strong 
interactions. The answer, at least for lattice QCD, seems to be 
positive IU9I1 though larger volumes and better statistics are 
still needed to further clarify the role of localization in QCD. 
In this final section we analyze why an Anderson transition is 
expected in the ILM and to what extent these arguments may 
also be valid in lattice QCD. 

As is known, unlike the zero temperature case where the 
decay is power-law, the fermionic zero modes in the field of 
an instanton at nonzero temperature (usually called a caloron) 
have an exponential tail e~ rT in the spatial directions and are 
oscillatory in the time direction B15I1 . This suggests that the 
overlap among different zero modes is essentially restricted 
to nearest neighbors in the spatial directions. However, in 
the time direction different zero modes strongly overlap due 
to the oscillatory character of the eigenmodes. This situa- 
tion strongly resembles a 4D disordered conductor in the tight 
binding approximation (only nearest neighbor hopping) with 
one dimension (time) much smaller than the rest so the sys- 
tem can be considered effectively three dimensional. It is well 
established that such a system may undergo an AT depend- 
ing on the disorder strength. In our case the role of disor- 
der is played by temperature since the wavefunction decays as 
~ e~ rT . From this discussion it is natural to find an Anderson 
transition in the ILM for a particular value of the temperature 
as first suggested in Ref . Jj] • 

The principal ingredient to reach the AT is thus an expo- 
nential decay of the eigenmodes explicitly depending on the 
temperature together with the possibility to tune the effective 
range of the exponential in order to reach the transition re- 
gion. Any theory with these features very likely will undergo 
an AT for some value of the parameters (see Ref. Il2tl for a 
more detailed discussion). It is therefore quite possible that 
the QCD vacuum has similar properties. Even if the objects 
responsible for localization are not the classical instantons this 
scenario could still be realized if the interactions behave in a 
similar way. We refer to [12] for a more detailed account of 
this interesting issue. 

In summary we have found that in the unquenched ILM the 
chiral restoration transition occurs around the same tempera- 
ture as the localization transition close to the origin. This in- 
dicates that both phenomena are closely related. From a phys- 
ical point of view this is not surprising since the appearance 
of the condensate is deeply linked II 1611 with the derealization 
of the zero modes in the instanton vacuum due to long range 
instanton induced interactions. It will be very interesting to re- 
peat these same studies in lattice simulations where the chiral 
transition also coincides with a deconfinement transition and 
additionally to see if the detailed structure of the topological 
objects affects the nature of the localization transition. 
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